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We study the persistent currents induced by the Aharonov-Bohm effect in a closed 
ring which either embeds or is directly side coupled to a quantum dot at Kondo reso- 
nance. We predict that in both cases, the persistent current is very sensitive to the ratio 
between the length of the ring and the size of the Kondo screening cloud which appears 
as a fundamental prediction of scaling theories of the Kondo effect. Persistent current 
measurements provide therefore an opportunity to detect this cloud which has so far never 
been observed experimentally. 



I. INTRODUCTION 



The Kondo effect has become one of the most studied paradigms in condensed matter theory in 
decadesJd The effect results from the interaction between conduction electrons and localized magnetic 
impurities. As a paradigm, the Kondo effect has played a considerable role in the development of 
renormalization group (RG) scaling ideas. The Kondo effect is associated with a large distance scale 
~ hvp/Tfc, where vp is the Fermi velocity and Tk, the Kondo temperature, which also defines the 
energy scale of the problem: Tk ~ ZJe"^/^, where D is the bandwidth and A the bare dimensionless 
Kondo coupling. The heuristic picture associated with this fundamental length scale is that a cloud 
of electrons with a size of order surrounds the impurity, forming a singlet with it. The remaining 
electrons outside the cloud do not feel the impurity spin any more but rather a scattering potential, 
caused by the complex formed between the cloud and the impurity, resulting in a 7r/2 phase shift at 
the Fermi energy. According to the renormalization group approach to the Kondo problem, all physical 
quantitieS|-siiould be expressed in terms of universal scaling functions of r/S^x, r being the distance to 
impurity .Bu Unfortunately, this large length scale of order 0.1 microns in metals has never been detected 
experimentally. 

Quite recently, due to some experimental breakthroughs, various aspects of the Kondo effect have been 
measured in a semiconductor quantum dot coupled via weak tunnel junctions to leads and capacitively 
to gates. By tuning the gate voltage, the number of electrons inside the quantum dot can be adjusted 
due to the Coulomb blockade. In particular, when the number of electrons is odd, the quantum dot 
can act as a localized spin S — 1/2 impurity. Nevertheless, unlike magnetic impurities in metals the 
physical parameters of the quantum dots can be varied continuously. When the tunneling matrix elements 
between the leads and the quantum dot are small, the system can be essentially described by a Kondo 
model. It has been predicted theoretically for a Kondo impurity embedded between two leads, that the 
transmission probability should reach one at low temperature T < Tk whereas it should be small at 
higher temperatur|fi. Itl Such a manifestation of the Kondo effect has been confirmed experimentally in 
the last few years. Ql 

The Kondo temperatures in these experiments are generally considerably smaller than I'^K and most 
importantly can be tuned via the gate voltage Vg. Therefore, the Kondo length scale is expected to be of 
order 1 micron or larger. This new experimental realization of the Kondo effect seems therefore to offer 
new opportunities to measure the screening cloud. In [ |l^, we have emphasized that the transmission 
probability of a quantum dot, in the Kondo regime, and embedded in a closed ring, may be sensitive to 
the length of the ring versus the Kondo screening length. The idea is simply that a finite size should 
suppress the Kondo effect even if the temperature is much less than Tk- Such an idea, central to 
the interpretation of the Kondo problem, has been for example used by Nozieres to interpret Wilson's 
numerical renormalization group approach of the Kondo problemEll. 

The purpose of this paper is to analyze in more detail the sensitivity of the transmission probability 
of the quantum dot to the Kondo length scale p. twp different geometries: one where the quantum dot is 
embedded in a closed mesoscopic cimjla.r ringE3~Ljl, and another one where the quantum dot is outside 
the ring and couples directly to it.E3'E3E3 We have shown schematically both devices in figure |l|. In 
both geometries, the crucial point is that the screening cloud is "trapped" in the ring and cannot escape 
into macroscopic leads. A natural way to measure the transmission probability through a quantum dot 



in a closed geometry would be by persistent current measurements. One might expect intuitively that 
the persistent current, as a function of the flux, $ penetrating the ring, will be much different when 
the screening cloud is small compared to the circumference, of the ring than when it is much larger. 
Such persistent Giir.rent experiments have been performed recently on micron sized rings not containing 
a quantum dotliHlla. We have shown that, when S,k/1 « 1, the persistent current is that of a perfect 
ring with no impurity (j cx evp/l) for the embedded quantum dot and is vanishing for the side coupled 
quantum dot. On the other hand, when S^k/^ >> 1, jl becomes much smaller in the first geometry, 
vanishing as a power of the bare Kondo coupling, whereas it converges toward that of a perfect ring in 
the second geometry. We always assume in this article that I ^ a, with a the lattice constant. In general, 
we expect to be a universal scaling function of £,k/1, in the usual scaling limit of the Kondo model 

(i.e. at small Kondo coupling and large ring size compared to the lattice constant). 

The plan of the article is as follows. In section 2, we first introduce the two tight binding models we are 
to study and develop a continuum limit analysis of the models in an open geometry which clarifies the 
difference between the two models. In section 3, we calculate the persistent currents for the embedded 
quantum dot in the Kondo regime in both limits £,k ^ I and ^ I, and discuss the extension of these 
results for asymmetric tunneling amplitude between the quantum dots and the wires. A special emphasis 
on the role played by particle hole symmetry is also placed. We also study in this section how bulk 
interactions modifies the conductance of the quantum dot and hence the persistent current. In section 
4, we present similar calculations for the side coupled quantum dot. Finally, in section 5, we discuss our 
results and the approximations we have used. Some technical details about our calculations have been 
relegated to three appendices. 




FIG. 1. Left: Embedded quantum dot. Right: side coupled quantum dot. Vg is the gate voltage and Vr/l, V 
control the tunneling amplitude between the quantum dot and the ring. 

II. FORMULATION OF THE MODELS 

The first model we are considering corresponds to a quantum dot embedded in a wire. In order to 
describe this situation, we begin with the simple tight-binding Anderson model: 

Hi=-t ^ (c]cj+i + h.c.) - t^(c]cj+i + h.c.) - t'[cli{c-i + ci) + h.c] + edC^Cd + Und]ndi. (2.1) 
i<-2 j>i 

Here sums over the electron spin index are implied. The quantum dot denoted by the subscript d simply 
corresponds to the site of the chain, t' , Cd and U are the tunneling matrix elements, gate voltage and 
charging energy of the dot, respectively. Uda is the electron number at the origin for spin a. 

The second model we want to consider corresponds to a situation where the quantum dot is outside of 
the ring. It has been referred to as the side coupled quantum dot. This can be modeled by the following 
tight-binding Anderson model: 

-t ^ + '^■c-) - t'lc^CQ + clcd] + tdc\cd + Und^nd[ , (2.2) 

i— — oo 

with similar notations and conventions. In this model the quantum dot couples directly to the site 0. 
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In order to distinguish both models, we wiU use in the foUowing the notations 1 for the embedded 
quantum dot (EQD) and 2 for the side coupled quantum dot (SCQD). These two models have already 
received-some attention in the past. The EQELhas been first studied in [ ^ and [ ||. In particular, Ng 
and LeeO have shown using Langreth formulallj that the conductance at zero temperature (and t' ^ 0) 
reads: 



Gi = — — sin (5, 

n. 

with 5 the phase shift which has been approximated by: 



(2.3) 



(2.4) 



and {ud) the average occupation number of the quantum dot. Nevertheless, we would like to stress that 
the phase shift according to the generalization of the Friedel sum rule by Langrethllj is proportional to 
the number of electrons displaced by the impurity "among which are included not only the d electrons, 
but also some of the conduction electrons" . The approximation of Eq. ( p.4[ ) should be valid only when 
the reservoir bandwidth is much larger than the bare Kondo coupling (proportional to t'"^) which implies 
that the contribution of th e cond uction electrons to the phase shift is negligible. We will come back to 
this point later (see section HID) with a special emphasis on the role played by particle-hole symmetry. 
Kang et aLEj have proved using similar reasoning that the conductance of the SCQD can be expressed 
as: 



Go 



2e2 



cos^ 5, 



(2.5) 



with similar notations. 

In both models, the Kondo regime corresponds to the limit: 



t' « -Cd^U 



Ed, 



(2.6) 



(e^ < 0) . In this limit the dot essentially is locked into the singly occupied state and thus corresponds to 
an S" = 1/2 spin impurity. In this regime, the EQD is predicted to exhibit perfe ct conductance according 
to Eq. (2^) whereas the SCQD should give zero conductance according to ( ^.5[ ). Note that at J = t' = 0, 
we obtain the opposite result. 

In the sequel, we want to recover these results from the continuum limit analysis of both models 
in the Kondo regime. In order to compare both models, it is instructive to write down explicitly the 
corresponding Kondo models. In this regime, the tight-binding Anderson model reads: 



i?i = ^ {c]cj+i +h.c.) + -i^(c]cj+i +h.c.) + J(c^ + c\)^{c-i + a) ■ S. 
j<-2 j>i 

Here the Kondo couphng is: 



(2.7) 



J = 2t'- 



1 



1 



-Ed U + ed 



(2. 



In general a potential scattering term is also induced, except in the case of the symmetric Anderson 
model, U — —2ed- However, provided it is small, it has little effect on the t ran smission. 

For the SCQD in the Kondo regime, the tight binding Anderson model (2.2) simply reduces to to the 
standard Kondo model: 



-t-oo 



H2 = -t^(cl+iQ + h.c.) + Jcl^co-S 



(2.9) 



At first sight, the two Kondo Hamiltonians (2.7) and ( |2.9| ) just differ in slight details. Nevertheless, 
when the Kondo coupling constant is switched off (J — 0), the first model reduces to two independent 
quantum wires and has therefore a transmittance T = whereas the second one has a transmittance 
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T — 1. Let us analyze the two models in the continuum limit. We can first linearize the spectrum around 
the Fermi points ±kp and introduce left (L) and right (R) moving chiral fields: 



ip{x) - e-''""='ipBXx) + e''""'i^L{x) , (2.10) 

■0 being the fermionic field associated with the operator c in the continuum limit. The next step is to 
introduce the even/odd basis defined as follow: 

V'e/oW = -^\Mx) ± i'Ri-x)]. (2.11) 

Observe that these two fields are left movers. It may be advantageous in some situations to formulate 
the Kondo Hamiltonians as boundary problems. 

For this purpose, we can fold the system by setting for a; > 0: 

^R.e/o{x) = l/je/oi-x). (2.12) 



All t he fields are now defined on the infinite half line x > 0. In both models, the Hamiltonians (2.7) and 



(2.9) take the form: 

H = Hf'^ +H°'""' + HZ^'', (2.13) 



with 



oo 



^cvcn/odd ^YL I dx 



2tt 



e/oi^)^^^'^L,e/o{x) - (ll^^j^^^^ ^{x)id^tp j^^^/ ^{x)) , (2.14) 







the Hamiltonians of free fermions. Notice that the Kondo interactions couple only to the even sector which 
is the main advantage of this even/odd basis. In order to derive the continuum limit for the interacting 
parts of both Hamiltonians, we need to pay attention to the boundary conditions at the origin before 
including the interactions. For the EQD, we had free boundary conditions on site 1 and —1 on the lattice 
implying therefore cq = 0. Conse quent ly, we have t/jjO^) ~ "0(0^) = and therefore tjjii{0~^) — —iPl{0^) 
and iPr{0~) = —-0^(0^). Using ( |2.12| ) and (2.11), we can infer that V'ii,e(0) = — V'L,e(0).r-Similarly, for 
the SCQD, we can prove that ?/'fi,e(0) — V'L,e(0) using the continuity of ipL/nix) at a; = 0.c2l 
The boundary conditions can be summarized as follow: 

EQD: ^j^,e(0) = -0L,e(O), (2.15) 
SCQD: Vi?,e(0)= ^L,e(0). (2.16) 

The interacting part of the EQD reads: 

HZ" = ^(^^1) + • ^(^(1) + V'(-i)) 

= g(Vl(l)e-''=- + 0l(-l)e^^--)| • S{M^)e^''' + M-l)e-"'n- (2-17) 
We then use the app roxim ations "06(1) = 0L,e(l) ~ ■0L,e(O) and ?/;e(— 1) = ^fl,e(l) ~ i'R.eiO) and the 



boundary condition (2.15). Finally, after these manipulations, the associated interacting parts of the 



Kondo Hamiltonians read: 



2 

where the bare Kondo couplings are defined as follow 



HZH = vfK^IM^ ■ Si^LAO), (2.18) 



^we have used the notations 0^ and to differentiate fermions to the right (+) or left (-) of the dot 
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4Jsin^fci? 2Jsmkp /r, in\ 

Ai = = , (2-19) 

TTVp TTt 

X2 = — = ./., ■ (2.20) 

Note that the derivation of ( p.20| ) comes simply from ^L,e(0) = ^(0)/\/2 for the SCQD. Therefore, both 
Kondo Hamiltonians have almost the same continuum limit, the main difference being in the boundary 
conditions at the origin. 

In both models the effective Kondo coupling renormalizes as: 

dX,/d\iiD = -X^ +0{X^) i = l,2. (2.21) 

D is the momentum space cutoff and effectively the bandwidth. In both models, we expect the effective 
dimensionless kondo couplings to renormalize to strong coupling namely Ai oo. On the lattice, this 
corresponds to an electron near the origin forming a singlet with the impurity spin. The remaining 
electrons can propagate freely except that they cannot go into the same local orbital as the screening 
electron since this would break the singlet and cost a large energy of order J. For the Hamiltonian H2 , 
the screening electron resides at site 0. Thus, at J — > 00, the other electrons cannot cross the origin; 
they are trapped on either the positive or negative axis, corresponding to perfect reflectance. Formally, 
the low energy effective Hamiltonian is simply: 

H2,iow = -tYl + ^-^-^ " + h.c). (2.22) 



Note that this effective Hamiltonian is similar to the embedded Kondo Hamiltonian (2.7) at J — 0. On 
the other hand, for the Hamiltonian Hi, the screening electron is in a symmetric orbital on sites 1 and 
— 1 {{ciy = {c\ +clj^)/\/2). The remaining electrons are now allowed to occupy sites 1 and —1 provided 

that they only go into the anti-symmetric state: cjj = (c°)^ = (c| — c^_i)/^/2. To obtain the low energy 
effective Hamiltonian in this case we must project the operators c±i into the anti-symmetric state: 

Pc±iP = ±-^. (2.23) 
The effective low energy Hamiltonian then becomes: 

Hi,iow =~tJ2 + ^-c-) - + '^■c-) " 7f ("''-2Ca + 4c2 + h.c). (2.24) 



j<-3 3>2 



A simple calculation shows (see Appendix B) that this effective scattering Hamiltonian has a transmission 
probability: 

T{k)^s\T?k. (2.25) 

Thus when kp — 7r/2, there is perfect transmission at the Fermi energy. Since the conductance is 
determined by T{kp), this model exhibits perfect conductance at half-filling. This corresponds to a 
transmission resonance for non-interacting electrons; it occurs at half-filling since the model has particle- 
hole symmetry in that case. A major difference between the strong coupling limits of Hi and H2 is the 
sensitivity to particle-hole symmetry breaking of Hi but not of H2 . While H2 exhibits zero conductance 
for any filling factor. Hi exhibits perfect conductance only at half-filling. In the limit of small bare 
Kondo coupling, the case relevant to experiments, particle-hole symmetry breaking effects arc small 
i n b oth Kondo models, of order J/D. Similar terms are indeed neglected in both formula (^.3[) and 



(2.5) when approximating d by ^{ud). We would like to emphasize that due to the fact that particle-hole 
symmetry breaking is an exactly marginal perturbation in renormalization group language, we can expect 
the conductance to be close to 1 for the Kondo Hamiltonian, Hi, at weak coupling and any filling factor. 
The speci al ro le played by particle-hole symmetry is not apparent in the approximation leading to eqs. 



(2.3) and (2.5). 
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In our analysis, we have also supposed free electrons in the wire. The situation changes considerably if 
we include Coulomb interactions and has been extensively studied in ref. [ J2^ In this case, particle-hole 
symmetry breaking becomes relevant for repulsive Coulomb interactions.E30 Nevertheless, the Kondo 
resonance can still be reached (and therefore the unitary limit in the embedded dot) but only for one 
special value of the gate voltage V* (and therefore e^), which would make the obser vatio n of the Kondo 



resonance a pretty delicate task. We will come back to this point latter in section [II E . The fact that 
the unitary limit has been successfully reached in [ || may be an indication that Coulomb interaction 
effects are very small in the leads. 

A similar analysis of the strong coupling fixed point can be reproduced using the field theory language. 
The strong coupling regime Ai — > cxd corresponds to a 7r/2 phase shift in the even channel. In terms of 
right and left movers, it changes the boundary conditions at the origin as follow: 



^-^,6(0) = -V'i?,e(0), Ai = , A2 ^ +00, 

V'L,e(0) = V'fl„e(0), A2 = , Ai ^ +00. 



(2.26) 
(2.27) 



Therefore, the boundary conditions in the even channels of both models get completely interchanged as 
it is schematically depicted in figure 0. 




ITI=0 



ITI=1 




X = 

2 

FIG. 2. Schematic flow for both models between the |r| — fixed point and the |r| = 1 fixed point. 



The strong coupling fixed point of one model corresponds to the weak coupling fixed point of the other 
and vice versa. In this analysis, we have neglected all marginal operators which correspond to particle- 
hole symmetry breaking operators. The transmission probability depends on the scattering phase shift 
of both even and odd sector as follow: 



(2.28) 



where tpo — 0, y A; and ipe varies from n/2 (Ai = 0, A2 — 00) to (Ai — +00, A2 =0). From this analysis, 
we recover the results obtained from the explicit analysis of the low energy effective lattice Hamiltonian 
at strong coupling. 



III. PERSISTENT CURRENTS FOR THE EMBEDDED QUANTUM DOT 



We consider in this section the persistent current in the transmission Kondo model (2.7). Thus we 
consider a closed ring of I sites (including the impurity site, 0) and apply a magnetic flux: 

$ = (c/e)a. (3.1) 

(We work in units where h = 1. e > is the absolute value of the electron charge.) In the Anderson 



model of Eq. (2.1), this corresponds to modifying the hopping term between sites j and j -I- 1 by a 
phase factor e'"^ such that J^j Q^j = ck- Which links carry the phase can always be changed by a gauge 
transformation, Cj — > e^'^^Cj . The total flux is, of course, gauge invariant. Choosing, for convenience, to 
put the phase factors on the 2 links connected to the impurity site, we obtain the Hamiltonian: 

1-2 

H = -tY,{c]cj+i + h.c.) ~ t'lclie^^'^/^ci-i + e'"^^ci) + h.c] + eJ^Cd + Und^m^. (3.2) 
In the Kondo limit this becomes: 
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1-2 

H = -t ^(cjcj+i + h.c.) + J(e'"/2c|„i + e-*"/2^I)|(e-*"/2^,_i + e^^/^ci) • S. (3.3) 

The zero temperature persistent current is determined from the flux dependence of the ground state 
energy: 

dEa dEo 

3 — — c— — = — e— — . (3.4) 
d<^ da ^ ' 

We expect that the persistent current is given by universal scahng functions of the variables a and 
l/£,K in the limit 1,(,k » a (where a is the lattice constant, which we have set to 1): 

lj = f{a,^K/l). (3.5) 

These universal functions depend on the parity of the electron number, N and also, of course, on wh ethe r 
we consider the EQD or the SCQD. This behavior is a rather immediate consequence of Eq. (^.4|), 
expressing the persistent current as the derivative of the groundstate energy with respect to the phase a 
and the applicability of usual continuum field theory scaling arguments. Since a is a phase, associated 
with the boundary conditions, we expect it to have zero anomalous dimension. The universal part of the 
ground state energy is expected to have the canonical dimension of 1/ (length), and thus the result follows 
from standard finite size scaling hypotheses. Ij may be expressed as a universal function of a and either 
the dimensionless renormalized Kondo coupling Xeff{l) or equivalently of the ratio S,k/1- A non-zero 
value of a/l leads to the appearance of various irrelevant operators in the effective Hamiltonian which 
produce corrections to Ij, down by powers oi a/l. This scaling form implies that j can be calculated using 
renormalization group improved perturbation theory, at large £,k/1 where Ae//(Z) << 1. It also implies 
that the current at small will be a universal characteristic of the strong coupling RG fixed point. 

We now want to calculate the persistent current in perturbation theory in J. When J = the impurity 
spin is decoupled from the rest of the chain which effectively consists of iV — 1 electrons on Z — 1 sites with 
open boundary conditions. The ground state at 0*'* order in J, consists of a product of the impurity state 
(which may be spin up or spin down) and the filled Fermi sea for the rest of the system. The 0*'' order 
ground state energy is independent of flux. The single particle energy levels have energy En — ~2t cos fc„, 
with kn = Trn/l, n = 1, 2, 3, . . . , Z — 1. The annihilation operator appearing in the interaction may be 
expanded in Fourier modes: 

X= (e'"/2ci +e-'"/2Q_i) = Wy^[e"/2sinfc + e-'"/2sinfc(Z- l)]cfe. (3.6) 

From this stage, we have to distinguish the cases N even and N odd. Here N is the total number of 
electrons, including the impurity spin. Let us first consider the case of an even number of electrons. 

A. Perturbation theory for N even 

For even N, the levels are doubly occupied for n < N/2 — 1. However, the state with k = ttN/2 
contains only one electron which may have spin up or down. Thus the ground state is 4-fold degenerate 
at J = 0. This degeneracy is split in first order in J. It is convenient to decompose x into the term 
involving k = kp = ■nN /2l^ which we label xoj plus all other terms, which we label x'. 

X, = y|sinfc^ [e-/2 - (-l)^/2e— /2j e,,. (3.7) 
Only xo contributes in first order perturbation theory. Keeping only this contribution, the Kondo inter- 



action in Eq. (3.3) reduces to: 



8t/ (7 

H„u = — sin^ikp) sin^ (5/2)4^ -c^^ • S, (3.8) 



where 
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a + kpl ^ a + ttN/2, 



(3.9) 



and hence: 



sin2(5/2= l-(-l) 



= 1 - i'-n^/2 



cos a 



/2. 



(3.10) 



Diagonalizing Hint in the degenerate subspace picks out the singlet state, with the first order ground 
state energy: 



Eq = constant — sin^ kp sin^ (5/2. 



Thus the persistent current, to 0{J) is given by: 

3eJ 



je{a) 



I 



sin kps'ma. 



(3.11) 



(3.12) 



We see from Eq. ( |3.9D that, at large /, the current oscillates in N with period 4 and that the current for 
N/2 odd is given by the current for N/2 even with a shift of a by tt: 



(3.13) 



Notice that Eq. ( 3.13| ) is indeed exact to all order of the perturbative theory ignoring corrections of 
0{a/l). Indeed, we can show that it is a property of the two different types of propagators involved in 
the calculation of the ground state at higher order of perturbation theory (see eqs. ( 3.18| ) and ([A^)), 
therefore a property of the ground state energy itself, and thus a property of the persistent current. 

To proceed to higher orders of perturbation in J, it is convenient to include the xo part of the Kondo 
interaction, in Eq. ( |3.8D in the unperturbed Hamiltonian, Hq- This Hamiltonian can be diagonalized 
exactly. The ground state has all levels doubly occupied or else empty except for the kp level. It has one 
electron in the kp-\eve\ which forms a singlet with the impurity spin: 



\s>= («L,tI^> -«L.iI^>)/^- 



(3.14) 



(Here the double arrows label the spin state of the impurity.) The remaining perturbation. Hint, has a 
term quadratic in x' and cross-terms containing one factor of xo and one factor of x'. The second order 
correction to the ground state energy can be written in the standard form: 



i / drT < s\H,nt{T)H,ntms > 



(3.15) 



The propagator for the field %' is, for r > 0: 



n=l 



[1 - (-l)"cosa] sin2(7rn/0' 



(3.16) 



where 



'2tcoa{nn/l) + 2t cos(7riV/2/). 



(3.17) 



At large r the sum is dominated by terms near the Fermi surface, n = N/2 — m, with m << N, so we 
may write: 



G(r) 



/-I 

.hpj: 

m—l 



4 • 2, 
y sm kp 



[l-(-l) 
1 



cos a 



, (r>0). 



(3.18) 



It can be easily seen that: 
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T < X'\r)x'{0) >= r < x'{r)x'\0) >= G(r) = e(r)G(|r| 



(3.19) 



where G(|r|) is given approximately b y Eg . ( 3.18 ), and e(r) — 1 if r > and e(T) = —1 if r < 0. 

We have calculated in the appendix Al je{ct) using (3.15) and the propagator (3.18) to second order 
in J, for N even and large I. We have ignored corrections of order 0{a/l), with a the lattice spacing. 
The result reads: 



j^{a) = [sma[X + X^\n{lc)] + (1/4 + In 2) sin 25] + 0{X^) {N even) 

where c is a constant of 0(1) which we have not determined. 
Let us now analyze the case N odd. 



(3.20) 



B. Perturbation theory for N odd 



In this case there is an exact two fold degeneracy of the ground state, which has S = 1/2. The 
unperturbed (J = 0) ground state has all free electron levels doubly occupied for k < kp = ■nN/2 and 
empty for k > kp. (We choose kp to lie exactly half-way between the highest filled and lowest empty 
levels.) Now first order perturbation theory vanishes since it necessarily gives states with a particle-hole 
excitation. Second order perturbation theory is now straightforward since the ground state is unique, 
once we specify a value of S^^^^i = 1/2 for example. The exact time-ordered Green's function is now, for 
r > 0: 



I 



(7V-l)/2 



[1- (-l)"cosa]sin2(7rn/0e«"^, (r > 0). (3.21) 



where: 

^„ = -2tcosT:n/l + 2tcosT:N/2l. (3.22) 
Again, at large r this sum is dominated by states near the Fermi surface so we can approximate: 



G(r) « y sin^ fcF ^ [1 - (-l)("-i)/2-mcosa]e-''^""('"+i/2)/' 

1 (-l)(^-i)/2cosal 



I 

2 ■ 2, 

y sm kp 



, (r>0). 



sinh{-KVpT/2l) cos]i{t:vpt/21) 
For T < 0, we sum over the unoccupied levels, with n — {N + l)/2 + m, m — 0, 1, 2, 



(3.23) 
In this case: 



T <xUr)xp{0)>^S^pGiT)^~^ ^ [l-(-l)«cosa]sin2(^n/0e«"^ (r < 0). (3.24) 



n=(N+l)/2 



At large |r|, this becomes: 



G(r) 



: sin^ fcf ^ [1 - cosale^-^^^^+i 



/2)/l 



m—O 



■ sin k p 



(-l)(^-i)/2 cosa 



sinh(7r?;FT/2Z) cosh(7rwi?T/2Z) 



, (r<0). 



(3.25) 



Thus G(t) is given by the same approximate expression, Eq. ( p.23[ ), for either sign of r. As usual, we 
may write: 



T < Xair)xl{0) >= -Sa.pG{-T). 
The second order term in the ground state energy is: 



(3.26) 
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^(2) _ 



2! 



a, 6 



dTG{T)[-G{-T)]. 



(3.27) 



If we use our large t expression, Eq. (3.23) for G(r), then we get an a-independent term whose integral 
diverges at r = 0. This is just an artifact of the large-r approximation to G(t) and in any event, this 
term doesn't contribute to the persistent current. The other term is a-dependent and finite: 



i^Q^"* — constant 



cos^ a 



3 J2 sin-* kp 
4/2 



dr- 



1 



cosh {ttvft/21) 



= constant 



cos^ a 



3 J2 sin"* kp 

IvpTT 



(3.28) 



Note that, in this case, the result does not depend on the parity of {N — l)/2 and is periodic in a with 
period tt. We have generalized this result to third order in J in appendix 
not been able to prove it to all orders of the perturbation theory. From t : 
then straightforward to infer the persistent current: 



A 2. Unfortunately, we have 



le ground state energy, it is 



jo{a) 



161 



sin 2a[X'^ + 2X^ \n{lc')] + 0{\'^) {N odd). 



(3.29) 



where c' is another undetermined constant. 

Several remarks about the persistent currents expressions (3.2C) and ( 3.29| ) are in order. First of all, 
we note a large parity effect- at small A the persistent current is much larger for even N than for odd N. 
Furthermore, as already noticed the periodicity is different in both cases. 

We have chosen to present the calculations using the lattice propagator and then consider the limit 
I ^ a, a the lattice spacing. We have also checked that wc ob tain the same results using directly the 
continuum limit representation of the Kondo Hamiltonian (3.3). The best strategy is first to unfold the 
chain with free boundary conditions in sites 1 and / — 1 in order to work with left movers only (the size 
of the chain becomes so far 21). In this representation, the linearized Hamiltonian (3.3) takes the form: 



H ^ ^ J dx iplid^Tpi 
-I 



Ai 

Vp — 

2 



-■S 

2 * 



^l(O)' 



(3.30) 



From this point, the perturbative calculations become very similar to what was presented in the text and 
similar results are obtained using the propagators involving left movers. 



C. Renormalization group arguments and the strong coupling limit 



In this section, we want to comment on the /-dependence of the persistent currents ( 3.29|) and ( 3.2C| ). 
In both cases, to the order in A that we have worked, the result has the form: 

&Vp 

je/o = —fe/o[Kff{l),a], (3.31) 

where Xeff{l) is the renormalized coupling constant at scale I: 

Ae//(0 = A + AMn/ + ... (3.32) 
If this result persists to all orders in perturbation theory, this would imply the scaling form for j we 



have anticipated on general arguments in Eq. (3_^) (which is valid in the limit of small bare coupling A). 
In particular, this implies that, at / << S,k, the perturbative result is valid. The finite size of the ring 
cuts off the infrared divergences of perturbation theory. On the other hand, when I is of order of or 
greater, perturbation theory breaks down. It is well-known from numerical renormalization group, Bethe 
ansatz and other calculations that one may think of the effective coupling constant as renormalizing to 
GO in the low energy effective Hamiltonian. Technically, this corresponds to a fixed point of the boundary 
renormalization group. Thus it is very reasonable to expect that the persistent current for I >> can 
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be obtained by simply taking the limit J oo. This is expected to be valid because the persistent current 
is determined only by the low energy properties of the Hamiltonian and these are correctly represented 
by the infinite coupling fixed point. 

It is a straightforward matter to calculate the persistent current for the tight-binding Hamiltonian of 



Eg. (2.7) at infinite J. In this limit we may simply use the low energy effective Hamiltonian of Eq. 
( 2.24 ). It is convenient to think of this as being a simple free electron model defined on {I — 2) sites 
labeled a, 2, 3, . . . , ^ — 2. (The impurity site is eliminated and the 2 neighboring sites, 1 and / — 1 have 
effectively collapsed to one site due to the projection onto the odd parity linear combination.) We can 
add a flux to the model by changing the phase of any hopping term. In particular, since the hopping 
term between sites I — 2 and a already has a reversed sign, it is convenient to change the phase of that 
term, resulting in the model: 

Hior, = -i 51 + '^■^•) - * 11(4^^+1 + ^-c-) - 4s(e'^"^"^cl2Ca + c^ca + h.c). (3.33) 

i<-3 j>2 v2 

We see that this is a standard potential scattering model with an effective flux of a + tt. The tt flux 
shift is a_cpnsequence of the projection onto the odd linear combination. As shown by Gogolin and 
Prokof'evo the persistent current at zero temperature and large I is completely determined by T{kp), 
the transmission probability at the Fermi surface, in such a potential scattering model (deflned at a = 0). 
From Eq. ( 2.25| ), at half-filling, kp — 7r/2, T{kF) = 1- For this value of T{kp), the persistent current is 



the same as for an ideal periodic ring. By summing over the energies of the filled Fermi sea,with momenta 
k = (2?™ -I- a)/l, n = 0, ±1, ±2, . . ., it can be shown that this gives, for M spinless fermions, an energy: 

Eo = -2t cos(H//) ^^'-^ ^ ^''^^ = constant -I- -^[d^ (M odd) (3.34) 
sm7r/( znl 

,, , ,,, sin 27r(M)// vp , ,9 n /„ 

= -2tcos([a- 7r]/0 ^ = constant -f [a - tt]^ (M even), (3.35) 

sinTT/t 27r( 

where \Q\ equals the principal part of the angle 6*, which lies between — tt and tt and jumps by — 27r at 
B = {2n+ l)7r. ie. 

[9]=e, {\0\<n) 

= e^2Ti, (7r<6l<37r), (3.36) 

etc. Thus EQ{a) has minima at a = 27rn, for M odd and at 27r(ri -t- 1/2) for M even. When we include 
electron spin, for non-interacting electrons, we find a 4-fold periodicity in the total number of electrons, 
N . If N is odd, then the numbers of electrons of opposite spin are necessarily of opposite parity, so the 
energy is the sum of the two terms in Eq. (3.34) and ( |3.35 ). On the other hand, if N is even, then the 



numbers of electrons of spin up and down are equal in the ground state, so the energy is twice Eq. (3.34) 
if N/2 is odd and twice Eq. ( plsj ) if N/2 is even 



To apply this free electron result to the large I (i.e. stron g coupling) limit of the Kondo model, we 
should take into account the shift, a — a -I- tt in Eq. ( ^.33 ), and also the fact that two electrons are 



removed from the low energy subspace corresponding to the impurity spin and th e sc reening electron 



These two effects actually cancel. Writing the result in terms of a defined in Eq. (3.9), we thus obtain 
the result ior I » ^k' 

. . . 2evF 1 /,r N 

je(a) = — a-TT, (iV even) 

TTl 

jo{a) = -^i[a] + [a~n]) (iV odd). (3.37) 

TTl 

The current is paramagnetic (i.e. positive for small positive $, corresponding to a current which 
increases the flux through the ring) for TV odd or for N and N/2 even. It is diamagnetic for N even and 
N/2 odd. This is exactly the same situation which we found at I « (^k in eqs: ( |3.20 ) and (3.29). It 



seems plausible that this feature is independent of S,k/1- We have plotted the persistent currents for N/2 
even and for N odd in Figures || and |^. 
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« / 7T 

FIG. 3. Persistent current of the EQD versus ajix for N = for ^,k/1 ~ 50 (solid line) and for ^,k/1 <C 1 
(dashed line). The persistent current for A'^ = 4p + 2 is obtained by a simple translation of tt. Note that 
Ae//(0 = l/ln(WO- 



0.5 



> 



-0.5 




O 



1 



3 



FIG. 4. Persistent current of the EQD versus q/tt for N odd and for for S,k /I ~ 50 (solid line) and for ^k/' ^ 1 
(dashed line), jo has been multiplied by x5 for visibility. 



D. Particle- hole symmetry breaking 

A crucial aspect of Kondo physics is the exact marginality of particle-hole (P-H) symmetry breaking. 
By exact marginality we mean that the P-H symmetry breaking coupling constants grow neither larger 
nor smaller under renormalization to all orders in perturbation theory. (By contrast the Kondo coupling 
itself is marginally relevant, since the corresponding RG /3-function vanishes in first order but is non-zero 
in second order.) 

The model is easiest to analyze when it has exact P-H symmetry: for instance the Anderson model of 
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Eq. (2) with = —U/2 at half-filling. Then various properties of the low energy physics follow exactly 
from symmetry. The average occupation number on every site < rii > has exactly the P-H symmetric 
value 1. This implies that the formation of the Kondo screening cloud has no effect whatsoever on the 
charge density. The peak in the single electron density of states (Kondo resonance) occurs exactly at the 
Fermi energy, E = 0. The phase shift at the Fermi surface is exactly 7r/2, as follows, for instance, from 
Langreth's application of the Friedel sum rule. Thus, the transmission probability at the Fermi surface 



is exactly 1 and the persistent current should have exactly the value for a perfect ring, Eq. (3.37) in the 
limit £_k/1 ^ 0. 

Once P-H symmetry is broken, all these properties change somewhat. However, the important result 
is that the changes to these properties are small, of order the dimensionless bare Kondo coupling. This is 
a non-trivial fact given that these are all low energy properties and that the effective Koado coupling at 
low energies is expected to renormalize to large values. A convenient way of seeing thisE3a is from taking 
the continuum limit of the Kondo model and then applying bosonization (either Abelian or non-Abelian) 
which leads immediately to spin-charge separation. One can then see that the Kondo interaction only 
involves the spin boson and the charge boson is apparently completely unaffected. If we start with an 
asymmetric Anderson model, Cd ^ —U/2, then the resulting large-?/ Kondo model contains a potential 
scattering term, of strength V, which is typically of the same order of magnitude as the Kondo coupling 
J. Under bosonization, this becomes oc Vdx4>c where 4>c is the charge boson. Since this is linear in it 
leaves the charge sector non-interacting and hence is strictly marginal. By factorizing the electron Green's 
function into charge and spin factors, it can be seen that this has the effect of changing the phase shift at 
the Fermi surface by an amount of order V/D. This is typically a small dimensionless number of order 
t'"^ /Ut ~ J/t Xq. There is a corresponding change in the electron density in the vicinity of the impurity. 
This, in fact follows from the Friedel sum rule as pointed out by Langreth. The phase shift at the Fermi 
surface is related to the total charge displaced near the impurity. In the Kondo limit < > is exactly 
1 and the small change in the phase shift is associated with a small change in the charge density at the 
nearby sites, i.e. in the charge density of the conduction electrons, rather than the "d-electrons" . Even 
if ed = —U/2 so that the Anderson model is nominally "symmetric", a P-H symmetry breaking electron 
density, different than 1, has a similar effect. In the Kondo limit, an effective potential scattering term is 
generated at second order in the bare coupling A. However, we again expect no large renormalization of 
this potential scattering term for the reason discussed above. In these cases the transmission probability 
at the Fermi energy, will be slightly less than 1 and the persistent current will be slightly modified from 



that of Eq. (3.37) even in the limit ^k/1 << 1- However, these modifications will be small, of order V/D. 
Thus, as a function of e^, for other parameters fixed, the transmission probability at zero temperature is 
expected to have a broad plateau where it is only slightly less than 1 for a range of of 0{U). Changing 
the filling factor only slightly modifies the shape of this plateau. 

On the other hand, things are very different if one considers P-H symmetry breaking in the Kondo limit 
with a bare Kondo coupling which is not small. In this case the potential scattering term V/D becomes 
of order 1 and the P-H symmetry breaking effects are large. This is true regardless of whether or not 
ed = —U/2. P-H symmetry breaking from the filling factor also has a large effect. Now the transmission 
probability is not, in general, close to 1 even when £,k/1 ~* 0. The modification of the charge density on 
the sites near the impurity is large, although (n^) itself remains at exactly 1 in the Kondo limit. In the 
extreme limit where the Kondo coupling goes to oo, it is ea sy to calculate the transmission probabili|tj{. 



and hence the persistent current, exactly. Using Eq. ( 2.25| ) and Gogolin and Proko'fev exact resulttll, 



the persistent current at arbitrary filling becomes in this limit: 



, . 2evF VTpsma / , r— \ r \ /t^t \ /r, 

je(Oi) — — = I arccos( yif cosq) — TTdjv 4p , [N even), (3.38) 

j,{a) = , yTpsina ^ arccos(yi> cos a) - n) , {N odd), (3.39) 

TTi V 1 — Tp cos^ a V ^ 



with Tp = T{kp) = sin {kp). When kp ^ 7r/2, these expressions differ considerably from (3.37) 



E. Inclusion of interactions in the ring 

Electron-electron interactions-jn_the ring have a dramatic effect on the conductance of the dot, and 
hence on the persistent current. c3'cll This can be seen by considering the strong coupling effective Hamil- 
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tonian of Eq. (2.24) for the EQD case. Now consider the effect of P-H symmetry breaking. For instance, 
if we choose half-fiUing, then we could include P-H symmetry breaking by a local potential scattering 
term at the origin: 

HlJon, ^ i?l,io„ + VcUa. (3.40) 

The corresponding continuum limit Hamiltonian is just the free Dirac fermion Hamiltonian defined on 
—00 < X < +00, with a back scattering term, V[iPI^{0)^Pr{0) + h.c] added. If V is small, this has little 
effect, reducing the conductance by an amount of 0{V'^). On the other hand, if we include electron- 
electron interactions in the rina^the potential scattering has a much larger effect. We can simply use 
the analysis of Kane and Fisheio (see also Ref . [ Q ) . Upon bosonization, and introduction of spin and 
charge bosons, we parameterize the strength of the (screened, short range) interactions by a parameter, 
gp related to the compactification radius of the charge boson by gp — I/ttR^. For zero interactions, 
gp = 2. Repulsive interactions decrease gp. (We assume that the interactions preserve the full spin 
rotation symmetry so that the corresponding parameter for the spin boson is unchanged from its non- 
interacting value, ga = I/ttR^ = 2.) The renormalization group scaling dimension of this backscattering 
term is: 

X = .gp/4 -I- 1/2. (3.41) 

This boundary interaction is exactly marginal for zero bulk interactions, relevant for repulsive interactions 
and irrelevant for attractive interactions. Thus, with repulsive bulk interactions, this back scattering 
parameter, V, renormalizes to large values at low energy scales corresponding to perfect reflection of spin 
and charge at the quantum dot. This implies that jl — > for large I. The consisiency of this assertion 
can be checked by considering the stability of the perfectly reflecting fixed pointEj At this fixed point, 
the ring is effectively severed at the quantum dot so that open boundary conditions can be applied to 
the fermion fields to the left and right of the dot. The leading irrelevant operator at this fixed point 
corresponds to a weak hopping process between the two ends of the severed chain, on cither side of the 
quantum dot. This operator has scaling dimension: 

x^l/gp + 1/2, (3.42) 

which is > 1 in the repulsive case. Hence this hopping process is irrelevant, for repulsive interactions, 
demonstrating the consistency of the assertion about the RG flow. We may think of this stable fixed point, 
where the Kondo coupling and the potential scattering, have both renormalized to oo, as corresponding 
to an electron in the even orbital at the original forming a spin singlet with the quantum dot with the 
odd orbital cither empty or doubly occupied depending on the sign of V. From this RG analysis-of the 
perfectly reflecting fixed point we can deduce the scaling of the persistent current with Icngth.EiJ The 
weak hopping amplitude scales with length as 

ie//(0«Zi/2-i/^''. (3.43) 

The conductance, being second order in ie//, scales with length as On the other hand, jl scales 

as trrff, going as the square root of the conductance, or the transmission probability at the Fermi surface, 
tJB . Thus, we expectEl: 

j (X Z'(i/fp+i/2)ginQ,. (3.44) 
This scales to zero faster than l/l for repulsive interactions. It is also interesting to consider the behavior 



for small V, close to the resonance. It follows from Eq. (3.41) that the effective back-scattering potential 
at length I is given by: 

Veffil) OC ?l/2-Sp/4^ (3 45) 

Thus the persistent current should depend on V and / as: 

jl:^ f{Vl^'^-^'/'^,a). (3.46) 

(This result only holds in the case << 1-) In particular, this implies that the width of the peak in 
j, as a function of V , scales as ?~i/2+£'p/4. Thus, for weakly repulsive interactions, gp < 2, j will have a 
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broad maximum as a function of gate voltage for intermediate I but this maximum will sharpen up with 
increasing ring length. The current at the maximum should be that of an ideal ring with no quantum 
dot. 

We note that, provided gp is not too small, and assuming symmetric leads, there is only one relevant 
operator at the resonant fixed point. Thus, tuning one parameter, such as the gate voltage, should be 
sufficient to pass through the resonance. Strict P-H symmetry is not necessary to get perfect conductance; 
even away from half-filling there should be one value of gate voltage where perfect conductance occurs. 
With asymmetric leads jl will scale to with I for all gate voltages. 

Further insight can be obtained by considering the weak Kondo coupling fixed point, with bulk inter- 
actions. Now the Kondo interaction breaks up into two parts which have different RG scaling dimension: 
the term that reflects an electron coming from the left or right and the term that transmits an electron. 
The folding transformation that we mentioned in Sec. II is not useful anymore because it would turn 
short range bulk interactions into inflnite range interactions. We use subscripts -I- and — to label the 
fermion fields to to the right or left of the dot. The open boundary conditions, at zero Kondo coupling, 
imply: 

V;i±(0) = -Vfl,±(0). (3.47) 
Thus the two parts of the Kondo interaction are: 

i^^„^ = VfX++ J2 ^i'L±-S + i^F A+- K+ I^L- + 4- fV-L + l ' S. (3.48) 



Initially 



2Jsin/ci? 

A++ = A+_ = A (3.49) 



However, while A++ remains marginal, even in the presence of bulk interactions, A-| does not. The 

marginality of A++ follows since the interaction can be written entirely in terms of the spin bosons (on 
the two sides of the dot) whose compactification radius (or interaction parameter, g^) remains unchanged 

by the SU (2) invariant bulk interactions. On the other hand, the dimension of A-| changes with bulk 

interactions. The corresponding operator has the same scaling dimension as a standard hopping operator 



between two severed chains, Eq. ( 3.42 ). Note that the spin-charge separation of the Kondo problem 
breaks down once bulk interactions are included. This results from the fact that we can't transform 

to even and odd channels when bulk interactions are present. A_| becomes irrelevant for repulsive 

interactions. We expect a flow from the A = fixed point either to the resonant fixed point if the 
backscattering term is tuned to 0, or otherwise to the charge and spin reflecting flxed point. The length 
scale at which the crossover to these other fixed points occurs is again given by « e~^/'^++ as before. 
This length scale should be the appropriate one to render the scaling variable Vl^^^^^p^^ , which controls 
the crossover between resonant and charge/spin reflecting fixed points, dimensionless. Thus we may write 
a general scaling form for the persistent current in the presence of bulk interactions: 

Jl = fi^K/l, {VlD){lliKfl''-3^l\ a). (3.50) 

Note in particular, that the width as a function of V is given by D(£^k /lY^'^^^''^^ ■ 



F. Asymmetric tunneling amplitudes 

In this section, we would like to study how the persistent currents are modified if the tunneling ampli- 
tude between the wire and the dot are non-symmetric. The tight-binding Hamiltonian (3.2) is replaced 
by: 



1-2 



H = -tY,{c]cj+i + h.c.) - [4(iLe^*"/'Q_i + tRe'^/^ci) + h.c] + eaclcd + Und^nai, (3.51) 
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where we have introduced the left and right tunneling amplitudes tL^tn. We have assumed no electronic 
interactions inside the ring. We refer to section [II E for a discussion of this case, the idea b eing that one 
more relevant operator needs to be adde d to take into account the channel asymmetry in (B.48). 
In the Kondo limit, the Hamiltonian ( |3.51 ) reads: 



1-2 



2 J 



tL 



tR 



tL 



-ia/2 I 



tR 



ia/2 
e ' Cl 



(3.52) 



where the Kondo coupling constant has been defined by J = {tj^ + tj^){^^ + jj^:^)- From ( 3.52 ), it is 
straightforward to reproduce our previous perturbative calculations. The persistent curr ent ( |3.20| ) and 
( 3.29| ) simply become: 



jeia) 
jo{a) 



Ksin<5[A + A2ln(/c)] + (l/4 + ln2)K2A2sin2<5] + 0(A^ k^) (A^ even). 



^'^^F('_^2 gjj^ 2a[A2 + 2X^ ln(/c')] + 0(A^ k'^) {N odd). 
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where we have introduced the ratio: 



'2tRtL 

j-2 I j-2 ■ 
+ ''fl 



(3.53) 
(3.54) 

(3.55) 



Notice first that the even and odd persistent currents are affected in a different manner. More importantly, 
the infrared logarithm divergences are only renormalizing the Kondo coupling as it should be. The 
persistent current is now a function of ci, Ae//,K. The universal scaling form for j (Eq. ( ^.31 )) is 
modified to: 



je/o{a,K,X,l) = ^Y~9e/o{£,K/l,a,K) 



(3.56) 



where the g^'s are universal scaling functions. Our renormalization group arguments prove that the results 
obtained using perturbation theory are valid only when S,k ^ I- 

In order to analyse what happens in the limit ^ Z, let us first observe that the asymmetric (infinite 
length) lattice Hamiltonian: 



Hasyrn = ^t ^ (cj^J + l + ^-C-) - t^iCjCj+i + h.C.) 
J<-2 3>1 



-2J 



tL 



tR 



■ s 



tL 



C-l 



tR 



V^l+^R 



Cl 



(3.57) 



has exactly the same continuuum limit (given by Eqs (^.13 , 2.14 , 2.18| )) as the symmetric Kondo Hamil- 



tonian defined in Eq. (2/7) (with J ^ J), provided we generalize the even-odd basis in Eq. ( 2.11 
to: 



tR 



tL 



^Lix) + 



tL 



tR 



i'R{-x), 



(3.58) 
(3.59) 



The generalized odd sector decouples from the generalized even sector which contains the Kondo in- 
teraction. We can therefore follow the same reasoning as in section 2. The effective Kondo coupling 
Ai = ^'^^'"^''-^ renormalizes as in Eq. (2.21). This effective dimensionless Kondo coupling is driven to 

strong coupling namely Ai — > cx). On the lattice, it means that a screening electron forms a singlet with 
the impurity in the orbital defined by: 
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= ( Jl.2 ^i + (3.60) 
In order to obtain the low energy effective Hamiltonian, we project on the orthogonal state Co'- 

Pc,P 



Fc_iF=-=^c„. (3.61) 



The low effective Hamiltonian then becomes: 



(3.62) 

It is straightforward to show that this Hamiltonian has the transmission probability: 

T{k)^{^%\\inH. (3.63) 



^R + ^R 



According to section HID, if t^jtn are small enough meaning the bare Kondo coupling small enough, 
particle-hole symmetry breaking effects are negligible. We can thus expect the transmission probability 

at the Fermi surface to be given by T{ep) — ^ . A similar result has been obtained in [^j^- In 

order to calculate the persisitent current, we can directly apply the resu lt of Gogolin and Prokof'ev.Ell 
The persistent current in the limit ^ / is therefore given by Eq. ( 3.39| ) with Tp = T{tF)- 



IV. PERSISTENT CURRENTS FOR THE SIDE COUPLED QUANTUM DOT 

We now want to consider the persistent current of a closed ring with Aharonov-Bohm flux and a side 
coupled quantum dot. As in the previous section, we suppose the ring has I sites. The Hamiltonian reads: 

i-i 

H = -tY,{c\+^c^ + h.c.) + ^ 4^ • ^^co. (4.1) 

i=0 

There is no explicit flux dependence in the above Hamiltonian since we can choose to gauge it away. 
Obviously, the flux dependence is encoded into the non trivial boundary condition Cj+; = e'"cj and in 
particular q = e*"co, with 4> — (cj e)a. Due to this a-dependent boundary condition, the possible values 
of the momentum are 

kn = 2TTn/l + a/l, 7i = -(/ - 2)/2, //2 (1 even). 



When J = 0, the persistent current is the one of free fermions given in Eq. (3.37) except that N has to 
be changed to N — 1. Note that the persistent current is an odd function of a and is 27r periodic for TV 
odd and tt periodic for N even. We will prove in the following that these relations persist to all order of 
the perturbation scheme in the Kondo coupling J. Let us first distinguish the two cases N even and TV 
odd. 



A. N even 

We first assume that a € [0,7r]. The ground state is well defined and has all levels empty or doubly 
occupied except the Fermi level at fciT- defined by 



17 



kp 



2tiN 

41 
2ttN 

41 



a 
J' 



N/2 even 
— a) 



I 



N/2 odd 



(4.2) 



which contains one electron forming a singlet with the impurity. Note that kp — kF{a). The interaction 
is the standard Kondo interaction between the quantum dot and the site of the ring. There is a 
contribution at first order which does not depend on a: 



p(i) _ 3J 



(4.3) 



The second order contributions can be obtained according to Eq. (3.15). Following the procedure 
developed in appendix A, we first expand ip{0) in Fourier modes and separate in the time ordered Green 
function the kp term (which we call -^o) from the rest of the series (which we call 'ijj'). Due to a dependence 
of kp, the left and right branches gives different contributions to the propagator. 

Let us first suppose that N/2 is even. We will show how the case N/2 odd is related to the case N/2 
even later. We have depicted a schematic level diagram in figure I for a e [0, tt] and N = 12. The lowest 
energy level has momentum k = a/l corresponding to n = 0. For N/2 even, the Fermi level is reached by 
one electron with spin up or down belonging to the left branch (corresponding to n = —3 on the figure). 

e(k) 







P 3 


-3 
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^ 1 




2 



kr^ 

P P 

FIG. 5. Level diagram for = 12 and a G [0,7r]. The filled circles depict filled levels with two electrons 
whereas the empty circles are for empty levels. The electron at the Fermi level (n=-3) has a spin I or J,. 

From this diagram, it is straightforward to compute the propagator for the field ip' which are taken at 
site 0: 



with 



G{t) = j e(^^-^-)^ ifr> 



(4.4) 



|fc|<feF 

1/2 



(4.5) 



where 



-2i cos[27r(iV/4 - to)// - a/l] + 2t cos[2tt N / 41 - a/l], 
-2t cos[27r(iV/4 -■m)/l + a/l] + 2t cos[27riV/4/ -a/l]. 



(4.6) 
(4.7) 



We assume that at large |t|, the sum is dominated by excitations near the Fermi surface. We can therefore 
use the approximation « —2TTmvp/l and « —2TTVp{m — a/vr), such that the Green functions can 
be written as: 
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Gi(T)ifr>0. (4.8) 



Similar calculations for t < leads to: 



G{t) = -j e(^--^-)^ ifT>0 

|fc|>feF 

' V 1 - e— r^^ 1 - e— r^^ 



-G2(r) if r < . (4.9) 



The first part of the propagator comes from the left excitations whereas the second part comes from the 
right excitations (see figure |^). 

We now want to emphasize that the only a dependence of the ground state energy comes from this 
propagator. Thus, properties of this propagator extend to properties of the ground state and therefore 
of the persistent current. First, notice that from time reversal symmetry it results that the ground 
state energy is an even function of a and therefore the persistent current an odd function of a. It 
is straightforward to show that the case N/2 odd is obtained from the case N/2 even by changing 



a ^ TT — a in (4.8) and (4.9). Note that it also a property of the Fermi momenta defined in (|4.2|). Using 



this property and time reversal invariance, we can prove that E^^_^_2{a) = {tt + a). We can go one 



step further. Indeed, it is worth noticing that under a ^ a + n the propagators (4.8) and (4.9) exchange 
with each other: Gi(t) — > G2(— r), G2(— r) Gi(r). Gi is associated with particle excitations and 
G2 wth hole excitations. When calculating the ground state energy in perturbation theory, each order 
of the perturbation theory appears as a time integral of a function which should be symmetric under 
interchanging Gi(r) with G2(— t) due to particle-hole symmetry. This proves the tt periodicity of (a) 
and also that: 



EilUa)=Eil\a). (4.10) 



We also need the correlator for ipQ-. 



T(cr|V'o,-y(T)V'0,i/(0)|e) = -jiSafSei, - 9{-T)Sy„Secr], 

T{a\iPoAr)^ljiO)\^) = \[-SajS,, + 0{t)S^J,^], (4.11) 

where |e) indicates an electron with spin e lying at the Fermi level. We have performed the calculations 
in the basis defined by: 

\A},\B) = I t,^),U,^), 

for convenience. The double arrows denote the spin state of the impurity. There are three contributions 
to the ground state energy at second order. The first contribution comes from the correlators involving 
only tp' and reads: 



00 

2; 



Attvf J ^2 2 





du G^{u)G2{-u)Tr -- [{A\S\u)S\Q)\B) + {A\S\Q)S'^{-u)\B)], (4.12) 



where we have made the change of variable r — 2ttvpu/1. We need to add also the cross terms (containing 
two "0' and two ipo). The full result for 

^-IJ dT{A\H.M{r)H,r.t {0)\B), (4.13) 

reads: 
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du [Gi{u)G2{-u) + -Gi[u)) Tr 



[{A\S°'{u)S\{))\B) + {A\S\Q)S''{-u)\B)] 



du (g^{u){A\A 



2 2 



cfc,5'^(u)5''(0)|i?)-G2(-^)(A|4.^ 



2 2 



+ j I du { G2i-u){A\cl ( ) Cfc,5''(u)5''(0)|S) - G^{u){A\cl ( ) c,,5''(0)5n-«)|i?) 



2 2 



where we have also used r = 2ttvpu/1. Using the identity: 



T{A\S-{u)S''mB) = i^'^'-^AB + '-eiuje-'^iAlS-lB), 
and gathering aU contributions, we finaUy obtain the second order correction to the ground state energy: 



du 



^G^iu)G2i~u) + ^^{G^iu) + G2i-u)) 



(4.15) 



Note that this expression is symmetric under Gi(u) <-> G2(— w) as it s hould be. It is also worth noticing 
that had we just focussed on the leading logarithm divergences in ( 4.15| ), we would find £^'■^^(0;) w 
~T7™7T^ ~ _ 3m^FA_in(/c) ^j^j^^j^ would rcnormalizc perfectly the first order contribution to the ground 
state energy given in ( |4.3| ) (it also provides a non trivial check of the calculations). Neve rthel ess, in order 
to calculate the persistent current, we are only interested in the a dependent terms in ( 4.15 ): 



S(2)(a) 



du 



3 [e^"^"^""/'^' + g-"(i+"/'^)] 
4 



9 



(1 



(4.16) 



This expression contains some UV divergences which are removed when we consider the persistent current 
definedby (r)(2) =-e£i?(2)(a): 



ij )l («) = 



du 



1 



(1 



1 



(4.17) 



The subscript 1 is for latter convenience. This expression is completely antisymmetric under a —^ it — a as 
it should be (since Eo{a) = Eo{7r — a)). In order to check explicitly the renormalization of the persistent 
current with the Kondo coupling, we would need to go the third order in perturbation. We have not 
performed the third order calculation but we expect the dimensionless Kondo coupling constant A = 
to renormalize as Xeff{l) = A + In Z + ... as in the previous section (see also further). Note that the 
a dependence of kp gives corrections smaller by a factor of a/l. This expression contains some severe 
infrared divergences when a — > or a — *■ tt. For example when a —^ 0, the level at fc^^ and the next 
higher level ( corresponding to n = 3 in Fig. ^ become degenerate, which is not taken into account in 
our initial choice of the ground state. For small value of a, we can evaluate the integral approximately 

by 



9ewi.A2 



eff 



321 



(4.18) 



Physically, this corresponds to truncating the summation over intermediate states occurring in second 
order perturbation to the first exited state only which approaches the ground state as a ^ and therefore 
gives the most important contribution to the persistent current. For perturbation theory to make sense 
we are restricted to <C ij'^)''^\oi)) namely 



a » TrAe// 
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(4.19) 



20 



When this condition is not satisfied, this approach fails. Note that this result is expected. Indeed, the 
persistent current is discontinuous at a = for A = which precludes a naive perturbative analysis near 
this singular point which corresponds simply to a level crossing. To overcome this difficulty occurring 
for a <C TT or |a — 7r| ^ tt, we need to do perturbation theory around the correct ground state which is 
built in this case from the two levels close to the Fermi surface which are mixed by the perturbation. 
In the sequel we consider the case a <^ tt, which extends trivially to |q; — 7r| <C tt by symmetry around 
a — -K 12. These two levels have momenta fci/2 = "F^^ + a/L They correspond on the figure || to 
the levels labeled by n = — 3 and 7i = 3. In fact, the strategy we follow is analogous to a second order 
degenerate perturbation theory. When a <C tt, the ground state is built with all levels with |fc| < k\ full 
and one electron lying on one of the two almost degenerate levels forming a singlet with the impurity 
(defining therefore two possible states). 

The first order correction in the Kondo coupling J mixes these two states. It is straightforward to show 
that at this order the contribution of these two levels to the ground state energy is found by diagonalizing 
the 2x2 matrix 

(^12) - ' 3/ . %J , (4.20) 



4; ^2 4/ 



written in the basis: 



|t,o;4)-U,o;^) \oA-A)-\oA-A) ..... 
^/2 ' V2 ' ^ ^ ^ 

where | tiO;-IJ.) = c| || |o,t;4J-) = CjtI Here the double arrows denote the spin state of the 

impurity. 

The ot her l evels are giving second order contributions i n J in ( 4.20| ). We have explicitly shown in 



appendix |C l| that the Kondo coupling constant in Eq. ( 4.2(]| ) gets perfectly renormalized at second 

/~i T r\ r\ T Tn /-iItt Til Vino +-<-\ Vn^ T/~ml n r \ /- 1 Vvt r n-/] ^ ^ f \ I \2l-,-,r7 ^1 \ / 7 n-o r. ^ f 1\ / 1 ' I ^Vi i T-z-icn 1 1+ no -n c^ it-t^t-t ct n /t 



order, namely J/l has to be replaced by 7ri;_F(A + ln[k])// = irvF^^ef f{l) /l- This result is not surprising 



since we expect the persistent current to be a universal function of Xeff and a as in Eq. (3.5). 
Therefore the final result for the ground state energy can be cast in the form: 



= + Ef-'^ia) + 2i sin[^] sm[a/l] - ^ sin'H^] + 

where £;^'''=^(a) is the ground state energy of the free case for N even. This is easily deduced from Eq. 



(3.37) by replacing iV by (TV - 1): 



^j/'^ee^Qf) = constant + |^ (a^ + (q, _ Tr)^) . (4.23) 

For Ae// = 0, the free case is recovered as it should be. The singularity in a = is smoothed by the 
square root. The persistent current becomes so far in the large / limit: 

•p/ X -e^'F / 2a ^ . ,,, 1 sin[2Q!/;i , 

Jl(«) = —r- 1 + cos[a/Z] - - , \ , (4.24) 

2./sin2M] + llA^,,^ 



-evF ( 2a 

TT 
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Notice that for Ae// 7^ 0, the persistent current is continuous in and (0) = 0. 
can be reproduced near a — n. Indeed, by symmetry it is enough to replace in (4.25) a/ir 



The same analysis 
a/iT 



1 



and to change the overall sign. It is clear that the result (4.25) is valid only for small value of a <^ tt and 
for Ae// <C 1 where the ground state can be regarded as a superposition of states belonging to the two 
almost degenerate levels. For larger value ofa_^ Ae// but still small compared to tt, we may evaluate 
the validity of this approach by developing (4.25) in powers of Ae///a: 



321 ^ 



(4.26) 



In this regime, we recover the small a limit of (4.17) 

Ae// < a <C TT 



le a range at sma 



Thus Eqs ( [4.17 ) and ( [4.25[ ) agree in the limit 

1 ^eff- 



We want to emphasize 



Therefore, together, they cover all t : 
that these arguments are valid for small value of Ae// ^ 1 meaning I <^ where perturbation makes 

sense. 

We have plotted the persistent current from Eq. (4.17) and (4.25) (more exactly j^l/evp) in figure ^ 
for S,k/1 ~ 50 (Ae// ~ 0.25). We notice that the amphtude of the persistent current is already strongly 
renormalized for this value of Ae//. Moreover the singularities at a = titt have been completely wiped 
out by the Kondo coupling. 

In the opposite limit / 3> ^k, where perturbation breaks down, we may expect the persistent current to 
be given by t he J 00 limit according to standard renormalization group arguments (see the discussion 



in section III). In this limit, we have shown in section 2 that the transmission T tends toward 0. According 
to the argument given by Gogolin and Prokof'evEll, we thus expect: 

j^^O for (4.27) 

the picture being simply that the ring is "cut" into two pieces which prevents non zero persistent currents. 



B. N odd 



The ground state is defined by all states under the Fermi level {kp — n{N—l)/4:l) being doubly occupied 
and all states above the Fermi level being empty. The number of electrons inside the ring is even therefore 
the system, including the imputity, is in a doublet state. Consequently, there is no contribution to the 
ground state at first order in J. In the following, we assume N = Ap + 1 and < a < tt. We will show 
later how the case iV = 4p + 3 is related to that case N ^ Ap+l. We have depicted a schematic diagram 
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of the different levels for = 13 in figure]^. The first level has a momentum a /I corresponding to ?i = 
and the last filled level has momentum k = (27rp — a) /I (labeled by —3 in the figure Q). 

e(k) 









P- 3 














r 1 




2 



F F 

FIG. 7. Level diagram for A'' = 13 and a £ [0, tt]. The filled circles are for full levels whereas empty circles are 
for empty levels. The Fermi level kp lays half between the highest filled level corresponding to n = —3 and the 
lowest empty level corresponding to n = 4 

The second order correction to the ground state can be written as 

+ C30 



p(2) 



dr T{0\H,nt{T)H,nt{Om- 



(4.28) 



The propagator reads: 



(4.29) 



where 



27r a 
4„ ~ —2t cos[—n + —\ + 2t cos kp 



.2tt a 
~2tcos[—n + yj 



2t cos 



27r(iV- 1) 



(4.30) 



As for A'' even , at large r, the sum is dominated by terms near the Fermi surface corresponding to 



^ . ^ — m, with m <^ N. Due to the a dependence of the quantum numbers, there is an asymmetry 



between left and right movers. By help of the level diagram it can be easily shown that 

1/2 ^ 1/2 



-T{m — a/27r) 



1/2 




.... ' + j2^- 

m—l ni—O 

^^r(l-a/27r) -{vFra/l) 



-r (m+a/27r) 



l-e- 



I 



1-e- 



(4.31) 



I 



This expression is valid whatever the sign of r. Following similar arguments as for N even, properties 
of the propagator translate into properties of the ground state energy. We can show explicitly that 
Eo{a) — Eo{a + 2tt) and especially Eo{a, 4p + 3) = Eo{a + tt, 4p + 1). Therefore the case A^ = 4p + 3 is 
simply related to the case A^ = 4p + 1 by a translation of tt. The second order term in the ground state 
energy reads therefore: 



CXD 

a,b ^ ^ _-L 



G(r)hG(-r)] 



(4.32) 
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After plugging the time ordered Green function (4.31) in this expression we get, retaining only a depen- 
dent terms: 



3.P 



SI{2'kvf) 



du 



1 + coshu(l — a/ir) 
1 — cosh u 



(4.33) 



This result can be generalized to the third order in J . Following similar calculations presented in 
appendix A2, we can show that: 



+00 



16 (27rwF)^ 



4nn(Zc) / dvG{v)[-G{-v)]. 



(4.34) 



We define the Kondo dimcnsionlcss coupling constant A as J = irvpX. Gathering all terms, we find 

+ 00 







'^''''^-[A2 + 21n[Zc]A3] 



16/ 



du 



1 + coshu(l — a/n) 
cosh u — 1 



(4.35) 







The correction to the persistent current finally takes the form : 

u sinhw(l — a/n) 



16/ 



du 



cosh u — 1 



(4.36) 







where we have defined as usual Ae// = A + A^ ln(/c) the rcnormalized labeled Kondo coupling constant. 
Exactly as in the case N even, this expression diverges when a ^ 0, i.e. when the two levels close to the 
Fermi surface become degenerate. In figure 0, it corresponds to the levels labeled by n = —3 and n — 3. 
For small value of a, we can evaluate the integral approximately by 



(j°)(2)(a)«- 



16/ 



Therefore, perturbation theory applies when {j°)'i\a) <C (ji)'^'^^'^) namely 



V3_ 
2V8' 



(4.37) 



(4.38) 



For smaller value of a ^ tt the Kondo coupling mixes strongly the states belonging to the two levels 
close to the Fermi surface. We need therefore to perform our perturbative calculations around a new 
ground state built from these two almost degenerate levels which we denote by 1 and 2. They have energy 
ei/2 = — 2/; cos(fcF T respectively, where we have defined kp — ^liE-Jl^ \Ye can now repeat the same 
analysis as for N even. 

The ground state has all levels with |fc| < ki full and two electrons in the levels 1,2. We have four 
possible states with Stot = S^^^ = 1/2: 



T,i;^) + U,T;^)-2| T,T;^) 

x/6 



T,i;^>-U,T;^) 



V2 



(4.39) 



The states are defined according to \a,(3,'i\) — c| ^Cj^l | Ti7o;f|") = c[ J fl-), etc. The first order 
contributions in J mixes these states. The associated matrix taking into account the contribution of 
these two levels to the ground state energy at first order in J reads: 



t J 



H12) = 



( £1 + £2 

J76 



J 



4/ ^'=1 

4/ " 




41 4/ 







262 

ei + £2 / 



(4.40) 
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We first note that the fourth state does not mix at tliis order and the matrix is effectively 3x3. We 
want now to include the effects of the higher levels which are giving second order contributions in J 
in ( 4.40| ). We have shown in append ix p 2| using second order degenerate perturbation theory that the 
Kondo coupling constant in the Eq. ( [4.40D gets perfectly renormalized namely J /I has to be replaced by 
T:vp{\ + ln[Zc])/Z — irypXef f /I- 

Defining E' = E — ei ~ €2 — the ground state energy is obtained by finding the minimum root 

|fe|<fei 

of the cubic equation: 



{E'Y + 



I 



{E'Y - E' 



4P 



TTVpX 



I 



(£1-62)' = 



(4.41) 



After the change of variable E' — ^^j^{X — and using €2 — ei ~ 2vfC(/1, the equation can be cast in 
the reduced form: 



, 



where a = 2a /n. 

The smallest root is found to be 



(4.42) 



Xmin ^ 2 ^/Q cos[^ + ^] with 6* = arccos 



-R 



(4.43) 



where Q, R are defined by: 

Therefore, the ground state energy reads: 



R 



2^ 3 



^(0) ^ j^frec + ,2 - ei + ^X„„„(2aA, Ae//), 

where E^'^'^'^ is the ground state energy for J = 0. 

The persistent current can be then readily deduced as: 



(4.44) 



(4.45) 



(4.46) 



with a = a ii (N — l)/2 is even and a = a — tt ii {N — l)/2 is odd. Moreover 



V 2a e 2tt 2VQ . 9 2n 

X^i„ = — = cos — sm — 6* (a) 

™" 3VQ ^3 3 ^ 3 ^3 3 ^ ^ ^ 



(4.47) 



and 



9'(a) 



\ r29\2 l„2^ 



(4.48) 



Notice first t hat f or A = 0, the free case is rec overe d as it should be. To analyze the range of validity of 
the formula (4.46), it is interesting to expand (4.46) in Ae//. It is straightforward to show that: 



(4.49) 



which correspon ds th e small a limit o f the persi stent current calculated with standard perturbation 
theory (see Eg. f|4.36D ). The refore Eqs ( |4.36D and ( f4.46D together cover all a for small Ae//. From both 
expressions ( 4.36| ) arid ( |4.46|) , we have plotted the persistent current for {N — l)/2 even on figure ^ for 
(Ae// «0.25). 
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« / 7T 

FIG. 8. Persistent current for A'^ = 4p + 1 calculated for S^k/I ~ 50 (solid line) versus the J = case (dashed 
line). The A'^ = 4p + 3 case is obtained by a translation of tt of the horizontal axis. 

Notice that the corrections are weaker for N odd that for N even. 

In the opposite hmit I, perturbation theory does not apply. However we can use the same 

reasoning as for N even in order to prove that 

f~^0 for (4.50) 

therefore no persistent current should be observed. 

V. DISCUSSION AND CONCLUSION 

In the two previous sections, we have calculated the persistent current in the two limits £,k/1 3> 1 and 
^fc <C 1 for both the EQD and SCQD. Our results are summarized in figures ^. In the embedded 

case, our results are actually very much different that those obtained in [ Q by solving some approximate 
self-consistent equations. In particular, the current was predicted to be small for N odd for large I/^k 
but not for small 1/£,k, the opposite of our result. Furthermore the tt periodicity in a that we find for 
odd N was not obtained even in the limit of small 1/^k where our perturbative calculations are robust. 
This variational approach seems to give the inverse behavior. We have, obtained similar disagreement for 
the SCQD with ref. [ and [ 0. We refer to our previous papertHi for a discussion of the validity of 
their results. 

In this article, we have analyzed the validity and robustness of our results against various perturbations 
which may occur experimentally. Let us summarize our results and also discuss some other potential 
experimental limitations like disorder. 

We have first analyzed how the results are modified if we consider electrons interactions inside the 
ring. From the theoretical point of view, particle-hole symmetry breaking terms are relevant for repulsive 
interactions. As a consequence, the gate voltage controlling has to be tuned to resonance in order to 
reach a perfect transmittanceB3 Therefore, it may imply experimentally that it might be difficult to reach 
the unitary limit and thus to observe the "saw-tooth" form of the persistent current which is expected 
for < I- 

We have also treated the case of asymmetric tunneling amplitudes ^t^, between the EQD and the 
wires. This affects especially the strong coupling regime I 3> £,k because the transmission probability is 
no longer 1 and therefore the persistent current loses its "saw-tooth" shape. Such effect could disguise 
the cross-over between the regimes I <^ S^k to / ^ From the experimental point of view, it would be 
important to tune the voltage gates Vi and Vr (see figure 1^) associated with and as close as possible 
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to each other in order to observe such cross-over. Moreover, in the presence of repulsive interactions, the 
operator corresponding to the tunnehng asymmetry becomes relevant at resonance. It would therefore 
imply experimentally one more parameter to tune. 

We have supposed the wire to be relatively clean such that we can neglect the effect of non-magnetic 
impurities. One may wonder how our results are modified when we have both the quantum dot and just 
one non-magnetic impurity with a potential scattering V localized at a distance r from the quantum dot. 
In the regime r < ^ ^ ^/j-, we can easily generalize our perturbative calculations for both the EQD and 
SCQD. We expect the persistent current to be given by a scaling function of a,£^k/l,r/l, V in this limit: 

jl^gia,^k/l,r/l,V), (5.1) 

where g is a universal function depending on the system we consider (EQD versus SCQD) but also on the 
parity of the system. We have explicitly checked for the EQD that the first correction to the persistent 
current is 0{JV^) for TV even and 0{J^V^) for N odd. In the other limit, I > r ^ (^k, it is reasonable 
to expect that the persistent current can be obtained by simply taking the limit J ^ oo. We can then 
easily calculate the transmission probability of the associated non interacting model. At half filling, the 
transmission probability reduces to the one of the non-magnetic impurity: T = Vp/ (vp + V'^) (away from 
half f illing, the expression is more complicated and has some r dependence). Plugging this result into 



( 3.39 ), provides the expression for the persistent current. The outcome is a decrease in the amplitude of 
the persistent current but especially the loss of the "saw-tooth" form here too. Such a result might also 
disguise a clear cross-over from ^ I to ^ I V is strong enough. From the experimental point of 
view, it seems quite important that the wires should be relatively clean. The fact that the unitary limit 
has been reached in [ |^ may indicate indeed that the wires are clean enough such that these effects are 
negligible. 

We have assumed only one electronic channels but despite the narrowness of present quantum wires, 
several electronic channels can be activated. However, it is reasonable to expect that one channel will 
have a stronger tunneling amplitude to the dot than the others. The RG equations, to third order, for 
the multi-channel Kondo problem are: 

j 

We see that if all but one of the couplings, Ai, are small, while Ai is larger and posi tive , then we may 
approximate the equation for the small couplings by only the second term in Eq. ( |3.2[ ), keeping only 
the AiA^ term. This equation then predicts that all the small couplings shrink. Meanwhile, the larger 
coupling grows. Therefore, the impurity spin is screened by an electron from the most strongly coupled 
channel and the other channels decouple at low energies. Thus, ioi I >> , we expect the single channel 
result obtained for the embedded quantum dot to still apply. Yet, we expect the results for the side 
coupled quantum dot to be fairly different. Indeed, if just one channel is screened by the spin impurity, 
the other channels do not feel the effects of the quantum dot in the scaling limit. The larger the number 
of channels, the less the transmission is affected by the dot. Therefore, if there are several active channel 
in a side coupled dot experiment, the observation of the cross-over between large £,k/1 to small 
could be severely concealed. In this respect, the embedded quantum dot would be a better candidate to 
detect the screening cloud through persistent current measurements. 

In this paper we have shown that the persistent current in a mesoscopic ring coupled to a quantum 
dot is a highly sensitive function of the ratio of the ring circumference to the size of the Kondo screening 
cloud. On the other hand, electron interactions in the ring, asymmetric tunnelling amplitudes, disorder 
and the presence of several channels can all serve to mask this sensitivity. Thus, while the screening 
cloud size is an important length scale in quantum dot physics, its separation from various other effects 
is likely to be a challenging experimental problem. 
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APPENDIX A: PERTURBATION THEORY FOR THE EMBEDDED QUANTUM DOT 



1. even 



In this appendix, we calculate the ground state energy at second order in perturbation theory for N 
even. The correction to the ground state energy is given by: 



=-4 / drT{s\xHr)^-Sx{r)xHO)^-Sx{0)U 



2! 



(Al) 



where x has been defined in (B.6). The first contribution to , coming from the term involving 2 
factors of x' in Hint, reads: 



(A2) 



where G(t) has been defined in ( 3.18 ). To evaluate the contribution of the cross-terms between xo and 
x' we need: 



^ < 7lxS,eWxo.i. (0)1/3 > = jsin^kFil-cosa)[S^,Sp^~0{~T)d,Jp^] 
T < 7|Xo,i.(t)Xo,, (0)1/3 > = y sin^ kpil - cos a) [-(5^,(5^^ + e{T)S,^6f!^], 



(A3) 



where I7) labels the state with one electron of spin 7 at the Fermi level. The full result for the ground 
state energy to second order in J can be written: 



a.b 



4 



s > 



G^ir) + - sin^ kpil ~ cos <5)G'(T)e(T) 



+ — siv? fci?(l — cos a) < s|c^^ 



(T (T 

T'T 



Ck,T[S''iT)S\0)\s>G{T) 



(A4) 



Here we have reinserted the creation operator at the Fermi surface, Ckp (at r = 0) in order to facilitate 
comparison with the first order term. Now using the impurity spin Green's function: 



this can be simplified to: 



T < S"^(t)5''(0) ^5"^ + ^e(r)e°'"=S'^ 



(A5) 



E. 



(2) _ 



G^{t) + y sin^ kpil - cos a)G(r)e(T) 



■sin^ A;F(l-cosa) < s\c\^'^Ckp ■ S\s > G{t)€{t) 



(A6) 



We now wish to examine the behavior of these various terms at large I. In the limit I » u_f|t|, 

4 sin^ kp 



G{r) 



TTVpT 



(A7) 



To estimate the large -I behavior of J dTe{T)G{T) we may use this expression for tq < |r| < l/vp where 
vpTo is of order t he la ttice constant (1 in our units). At larger r, G(t) decays to zero exponentially, as 
we see from Eq. (3.18). Thus, 
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dTe{T)G{T) 



) sin kp 

TTVp 



[In Ici + cos (5 In 2], 



where ci is a constant of 0(1), independent of I at large I. Now consider: 



G(r) 



16 sin kp 



dr 



2 cos (5 



(girtiFlrl/i _ 2)2 (g27r„F|r|/; _ 2) (e'^"f' 1^1 /' + 1)2 



(A8) 



(A9) 



The integral of the first term is divergent at r = 0. This just reflects our use of the large |r| form of 
G(t). Using the exact expression would give an integral of 0{P), corresponding to a term of 0(1) in 
the ground state en ergy . This is independent of a however, and so does not contribute to the persistent 
current, using Eq. (|3.4[). Thus we have: 



G(t) 



16 sin kp 
I 



2 cos a 



In/ 



finite 



cos^ a 



TTVp 



(-1 + 2 In 2) 



constant. 



(AlO) 



Importantly, the \nl terms in the term proportional to trcr^fi'', in Eq. (A4), cancel. Now evaluating the 
various terms in Eq. (A6), we have: 



E, 



(2) 



J^sin kp 
2\It:vp 



[cos a X finite — 6 cos^ a] 



-(1 — cosQ:)32(ln/ci + ln2cosQ;) < s|cj, ^c^^ • S\s >\ + constant. 



(All) 



Combining this with the term of 0(J), gives: 



En — 



4/ 



[cosa[A + A^ln(/c)] + cos^ 5(1/4 + In 2)A2] + constant, 



(A12) 



where c is a dimensionless cons tant which we have not calculated explicitly, and A is the dimensionless 
Kondo coupling defined in Eq. ( 2.19|) . 



2. iV odd 

Let us write in this case an expression for the term in the ground state energy of third order in J: 



f'^"' - ^2 5: tr (^^^) r < S'^{n)S\r,)S'^{r,) > 



■ J dTldT2dT3G{T2 - ri)G(T3 - T2)G(ti - T3). (A13) 

The factor of 2 arises from the two ways of ordering the three vertices. Here /3 is the inverse temperature 
which must be taken to oo. The integrals run between ±/3/2. We now use: 

T < 5"(ti)^^(t2)5^(t3) ^e'^''^e(Ti,T2,T3), (A14) 



where e(ri, T2, t^) = 1 if ti > r2 > and is completely antisymmetric in its arguments. Using Eq. ( 3.23 ) 
for G(t), we now get four terms proportional to cos" a with n = 0, 1, 2 or 3. The n = term does 
not contribute to the persistent current so we ignore it. The n = 1 and n — 3 terms vanish because the 
integrands are antisymmetric. This leaves only the cos^ a term, which is: 

.3) ^ -3JWa , ,^ r e(ri,T2,0) 



En ~ 5 sin*' kp dTidT2 -, , ^ , ; ^ —. (A15) 

° 2/3 J ^00 cosh{TrTiVp/2l)cosh{TrT2Vp/2l)smh[Tr{T2-Ti)vp/2l] ^ ' 

This integral is actually divergent at ti = T2, signaling the breakdown of the large t approximation to 
G{t). Integrating over vp\ti — T2I < / gives a contribution to the energy proportional to In/. We may set 
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Ti equal to Ti inside the argument of the (non-singular) cosh function, so that the logarithmic term can 
be obtained as: 

^(3) SJ^cos^a . e, 1 f ^ eW2Z 

EX' w 71, sm^'fcF / dTi ^ / dr \\ , (A16) 

where the r integral should be taken over tq < |r| < l/vp, where tqVp is a distance of order a lattice 
spacing (one in our units). Thus we obtain: 

^(3) 24J^cos^Q;sin'^/cFln(Zc') . 

E\' = -. — i — + constant, A17) 

(vftt) / 

where c' is another constant of 0(1) which we have not determined. Combining this with the second 



order result, of Eq. (3.2S) we obtain: 



Eo = cos^ ['^^ + '^^^ ln(/c)] + constant. (A18) 



APPENDIX B: J ^ cx> TRANSMISSION PROBABILITY IN THE EMBEDDED QUANTUM 

DOT 

We can easily calculate the exact transmission coefficient for the model with a + vr = in the Eq. 



(3.33). We simply solve the lattice Schroedinger equation: 

Ec^, = -i(0j+i + 0,-i), (j = 3, 4, . . . Z - 3) 

E4>2^-{t/V2)(l)a-t(f>3 
E,Pl^2 = -t<pl-3 - {t/V2)(f>a 

E(ba^~{t/V2)i<j)2+ (1)1-2). (Bl) 

The phase shift results from the reduction of the hopping matrix element by a factor of 1/ V2 at the site 
a, again a consequence of the projection onto the odd state. We can write down the general solution of 



Eq. (Bl), by inspection: 

0, = Aeij)smk{\j\ - 1) + Ssinfc|j|, (|i| > 2), 

(pa = BV2smk, (B2) 



with e{j) = 1 if j > and —1 if j < 0. In Eq. (B2), we label sites to the left of the impurity by negative 
integers, so 4>j is the wave- function at site I + with j = —2,-3,.... The allowed values of k depend 
on I but we don't need to consider them explicitly since we just need the transmission probability. For 
suitably chosen A and B we can obtain a scattering solution of the form: 



(f>j = e"'^ + VRe'^'e-'''^ (j < -2) 

4>j^VTe''fe'^\ (i>2), (B3) 

with 

T{k) = sin^ k. (B4) 

APPENDIX C: PERTURBATION THEORY FOR THE SIDE COUPLED QUANTUM DOT 

1. Degenerate perturbation theory for N even 

When a <C TT, we need to perform our perturbative calculations around a different ground state which 
is composed of the two levels near k^ = ttN/AI (we have used k^ which is different from kp defined in 
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Eq. (4.2)). In this section, we want to include higher levels in Eq. ( 4.20| ) at second order in perturbation. 
We make the approximation that the levels 1 and 2 are almost degenerate with energy ei « 62 ~ cm- 
We proceed strictly speaking with a second order degenerate perturbation theory. In other words, when 
introducing the sum over intermediate states (with energy e^) in the standard second order perturbation 
formula, we are approximating efc — ei/2 as — cm- We expect this approximation to be valid as far as 
|ei/2 — Cm I ^ — ea/ | which implies a <C tt. Obviously A is supposed to be small compared to 1 for 
perturbation to make sense. There are no contribution at first order in A from the other levels. The first 
contribution comes at second order: 



dT{p\TH„,tiT)H,, 



mi) 



(Cl) 



where p,q = 1,2 labels the two different singlet states (see Eq. (4.21)). Note that in our approximation, 
we suppose that the states \p) and \q) have the same energy em which is also a priori required for 
this formula to make sense. In the mode expansion of ip, we separate as before ijji/2 (corresponding to 
momenta ki/2) from other modes (which we call ijj'). The propagator for the ip' field remains unchanged: 



(C2) 



We also need the other propagators: 

^(*7lV'fe,c('r)V'fe',<5(0)|j/?) = j[Skk,Sk'kjSe^Ssi3 - Q{-T)5kk'Sk,kjSe55ji3] 

^(«7l'0fe',5(T)V'fe,(:(O)|j/3) = ji-Skk.Sk'k^SeiSsfs + QirjSkk'Sk.kjSesSjfs] 



(C3) 
(C4) 



whith k,k' = ki,k2, and where the electronic state [17) corresponds to on electron with spin 7 with 
momentum ki = fci/2- 

It can then be shown easily that 



41'=-'^ I drY^\trY-<p\S^S'\q> 



a a 



+ ]Gir) {P\cl 



a a 

T' '2 



CkT{p\S'^{T)S\m- 



(C5) 



Notice that the log(Z) terms in the first line of Eq. 
large I is: 



[Cq) cancel. We find that the leading behavior at 



pq 



3J2 



In Ic — 



In^c, 



(C6) 



which is independent of p, g and renormalizes perfectly the matrix elements of Eq. (4.20). We have 
therefore proven that in the range a ^ tt, the main effect of the higher levels is to renormalize the Kondo 
coupling, the infrared divergences being cut off by the size of the ring. 



2. Degenerate perturbation theory for N odd 

When Q ^ TT, the two levels near kp — ■k{N — 1)/4Z beco me almost degenerate (see Fig. l7|). We want 



now to include the effects of the other levels in the Eq. ( 4.40 ). They are contributing to the kinetic energy 
but are also giving second order contributions in J in ( [1.4C ) . We follow exactly the same scheme as for 



N even. Namely, we assume that the two levels 1 and 2 close to the Fermi surface are almost degenerate 
and have energy ti/2 = — 2icos[/cF]- By analogy with second order degenerate perturbation theory, we 
need to compute: 

E[f dT{z\THUr)H,ntmj)- (C7) 
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where 1 < i,j < 4 correspond to the four states defined in (4.3£) which are taken to have the same ground 
state energy ti ^ ep for ( |C7| ) to make sense. It seems more convenient to first compute the quantities 
A^^i = — i J dr {Sk\Hint{T)Hint{0)\Si) where the states \Sk) define the basis 

l^i) = I T,T;4) ; 1^2) = | T,i;^) ; l^a) = I i,T;^) ; 1^4) = I Ti,o;^) ; 1^5 =)\oAi-A) 

As usual, we separate in the Fourier decomposition of -0^ = e^'^^ipk the modes involving fci/2 = 

kp^a/l from other modes, which we label by V''- The propagator involving the prime fields only is easily 
performed: 

r(V.'t(.)0'(O)) = 2cosh(z;W0 



Let us define the electronic states \ri) and {rj\ as follow 



follow: 



(C8) 



where the momenta pi , p2 , 91 , 92 take the values ki or k2 and 7i , 72 , /3i , are the associated spins (the 
vacuum |0) means the Fermi sea i.e which has levels with |A;| < ki full). Let us also define the impurity 
state \ui) such that 

l^.> = h>(g)l".> 
The other propagators we need are the following: 



+Sk'p2^l^t2{^kq2^PiqiSel32^-tll3i Skqi5p^q^5el3^5^il32) " © ''')^fcfc' ^ei/^s; ] 

= y [-4'pi^i.7l('^fc'Zl'^P292<^£/3l^72/32 - '^fc«2^P2gi'5e/32'^72/3i) 

— '^fe'p2'^'^72('^feg2'^Pl9l'^e/92'^7l/3i ~ ^kqi^p^q^S^li^Sy^jj^) + © (r)(5fefc/ (5^^ (5s . ] 



(C9) 



(CIO) 



where A:, fc' = ki,k2- Using the expressions of various propagators, it can be shown after some algebra 
that 



A 



(2) 
SiSi 



2 



\Uj > 









'^P2 92 '^72/32 
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2 



^P2 91 ^J2f3l 



T'T 



/3i7i 

<5)9l(j2'^7l& 



I 



(Cll) 



0211 



'a" 
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/5i72/ 



/32 72 

r(u,|^''(r)5^(0)|i.,) 



We are only focusing on terms leading to infra red divergences which will renormalize the Kondo coupling 
constant. Notice that the first line of ( Cll ) do not give any logarithm contributions. The leading 
contribution at large I can be then computed straightforwardly in all cases. We can establish 
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(2) 
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(2) 
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(2) 
24 



2nvpl 



■ Inlc 



21 
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(2) 
13 



2t:vfI 



\nlc- 



21 



In Ic 



2lTVpl 



21 
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(2) 
34 



-A 



(2) 
25 



-A 



AlTVpl 



■ In^c 



TTVpX^ 



InZc . 



(C12) 
(C13) 
(C14) 
(C15) 
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All other matrix elements do not give an y lo garithm divergences. From these results, it is then easy to 
compute the matrix elements defined in (C7) and to show that the Kondo coupling constant in the Eq. 
(4.40) gets perfectly renormalized namely J /I has to be replaced by ttvf{\ + ln[Zc])/L 
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